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Abstract 



We explain why the maximal positive Lyapunov exponent of classical 
SU(iV) gauge theory coincides with (twice) the damping rate of a plasmon at 
rest in the leading order of thermal gauge theory. 
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I. INTRODUCTION 



Numerical studies of Hamiltonian SU(iV) lattice gauge theory in (3+1) dimensions have 



value of the largest positive Lyapunov exponent Ao has been obtained for SU(2) and SU(3) 
with the result |[g 



where E p is the average energy per plaquette, c 2 ~ 0.17 for SU(2), and c 3 ~ 0.10 for SU(3). 
For the SU(2) gauge theory the complete spectrum of Lyapunov exponents was obtained on 
small lattices |§. These calculations, which follow the evolution of a classical gauge field 
configuration in Minkowski space, also showed that the energy density distribution on the 
lattice rapidly approaches a thermal distribution [[| . This finding confirms the expectation 
of a finite growth rate of the coarse-grained entropy density of the gauge field, which follows 
from the observation that the sum over all positive Lyapunov exponents at fixed energy 
density grows like the volume jj. Hence, at any given level of coarse-graining, the classical 
gauge field "self-thermalizes" on a time scale of the order of the inverse Lyapunov exponent. 

In order to determine the value of the maximal Lyapunov exponent Ao, the evolution 
of the gauge field configurations must be followed over periods t 3> A^ 1 . The Lyapunov 
exponent is therefore effectively obtained for gauge fields that are members of a thermal 
ensemble, and we can identify the average energy per plaquette E p in (1) with that of a 
thermalized lattice. At high temperature the gauge field is a collection of weakly coupled 
harmonic oscillators, hence the average energy per independent degree of freedom of the 
classical gauge field is equal to the temperature T, yielding E p = |(iV 2 — 1)T for SU(iV). 
The factor | accounts for the restrictions imposed by Gauss' law. We can therefore rewrite 
the result (0) as 



shown that the gauge fields exhibit chaotic behavior in the classical limit ffl]. The numerical 



Ao = c N g 2 E 1 



(1) 




2 T (N = 2), 
2 T (N = 3). 



(2) 
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As already noted in [|J these values for Ao coincide, apart from a factor 2, with those of 

the damping rate of a thermal plasmon at rest, obtained by Braaten and Pisarski |J in the 

framework of thermal perturbation theory: 

N (O.UQg 2 T (iV = 2), 

7o « 6.635— g 2 T = (3) 
247r [0.26VT (N = 3). 

The goal of the present article is to establish this connection and to explain the origin of 

the factor Ao/70 = 2. 

We approach this goal in several steps. First we review the numerical "measurement" of 
the Lyapunov exponent in classical lattice gauge theory. We point out that the exponential 
growth rate of a small perturbation in the magnetic energy density used in those calculations 
is equal to twice of the growth rate of fluctuations in the elementary field variable, in the 
continuum limit the vector potential. This explains the factor 2 between Ao and 70 . 

In the next step we demonstrate that in classical calculations the linear perturbation 
propagation corresponding to the equations of motion of a chaotic dynamical system has in 
general a Fourier spectrum of imaginary frequencies. The Lyapunov exponent is equal to 
the magnitude of those imaginary frequencies. 

Then we argue that the chaotic dynamics of the classical system acts like a thermal 
ensemble averaging the perturbation propagation equation over stochastic frequencies. The 
square of these frequencies can either be positive or negative. In this case the damping rate 
and the plasma frequency of the classical elementary field fluctuations are related to the 
mean value and the width of the probability distribution of frequency squares. 

The final result of these considerations is that the Lyapunov exponent as defined in [1] 
measures twice the damping rate of classical gauge field fluctuations on the lattice. It is left 
to show that the quantum field theoretical calculation of the thermal damping rate at rest in 
hot perturbation theory in the leading 0(g 2 T) order survives in the classical (h — * 0) limit. 
We begin with the discussion of this point in order to establish connection with thermal 
quantum field theory. 
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II. COLLECTIVE PLASMA MODES 



We begin by briefly reviewing the derivation of the plasmon damping rate. Nonabelian 
gauge field fluctuations in a thermal background have been studied extensively in the frame- 
work of perturbation theory [|6|-|T0|j. The gauge field develops massive collective modes (plas- 
mons) with frequency u(k) > k due to interaction with "hard" thermal gauge bosons, i.e. 
excitations with energy of order T. The energy of a plasmon at rest is m g =w(0) = \\fNgT 
in SU(iV) gauge theory. For our purpose it is important that the dispersion relation u(k) can 
be obtained in the framework of semiclassical transport theory, where classical field fluctua- 
tions a M are coupled to the quantized thermal excitations of the gauge field The gauge 
invariant description of the collective modes requires the introduction of effective n-point 
vertices ||, which can be systematically derived from the effective action fT2| , [T3ll : 



AttlK) = -\m] J dh tr [r^^fi) (4) 

where n a = (l,n) is a null four-vector, and the integral is over all directions of the spatial 
unit vector n. D u stands for the gauge-covariant derivative. We have denoted the collective 
gauge potential and field strength f^ u by lower-case letters to indicate that these describe 
fluctuations around a thermal background. Note that £htl is a classical construction, 
with the sole exception that the plasmon rest mass m g depends on the energy distribution 
n(uj) = (e huJ / T — 1) _1 of quantized thermal excitations of the gauge field: 



2 2 AT f d 3 k , / •> \ Ng* 2 



Ng 2 — / f^n(w)(l + n(w)) = — ^-T 2 . (5) 
9 3 T J (2tt) 3 K JK K JJ 9 h K ' 

At leading order in g, (|5]) is evaluated for hard thermal quanta with u — \k\. 

Braaten and Pisarski || showed that the collective plasmon modes are unstable due 
to the effective interaction (|J). The plasmon damping rate 7(/c) is defined as imaginary 
part of the plasmon pole in the Feynman propagator corresponding to decaying plane wave 
solutions. The rate of instability for a plasmon at rest can be expressed as the imaginary 
part of the polarization function of the gauge field at the plasmon pole fT4|j : 



4 



7o = 7(0) = *n t (m 9 + «0, 0), (6) 

where the transverse polarization function *U t (uj,k) only depends on soft modes described 
by (0j). The plasmon rest mass exactly cancels from the expression (||) and the result ([3]) 
is a pure number multiplied by g 2 T, which is a classical inverse length scale. In fact, the 
calculation explicitly makes use of the classical limit of the Bose distribution, n(oS) — > T/hui, 
in the evaluation of the loop integral (see eq. (23) of ref. ||). 

' " , assuming 



Since the effective action can be derived from classical considerations |15[ 
a given spectrum of thermal excitations, it also applies to the collective excitations of the 
classical gauge field on a lattice. The sole modification is that the spectrum of thermal 
fluctuations is now given by the limit of the Bose distribution. Denoting the lattice spacing 
by a we find 

m 2 ^ 2 1 = J_ T 

9 3 y ^ uo 2 3vr y a V ; 

k 

in the weak-coupling, large volume limit. The plasmon mass (0) is a purely classical quantity 
of dimension (length) -2 not containing h, but it diverges in the continuum limit a — > 0. This 
is not surprising, since the lattice spacing serves as a cut-off that is required to regularize 
the ultraviolet divergences of the classical thermal gauge theory. The exponential growth 
rate of small classical field fluctuations is not affected by this divergence because it does not 
depend on the value of m g , as mentioned above. The result @ for the plasmon damping 
rate 70 remains valid if the correct plasmon mass m g in the effective action (|4]) is replaced 
by the value (|7|) for the classical gauge field defined on a lattice. 

More intuitively, the independence of 70 from the value of m g can be understood as 
follows. The cross section for scattering of a thermal gluon on a slow plasmon is: 

where /xd = V3m g is the inverse Debye color screening length. The scattering rate v is 
obtained by multiplying with the gluon density in the initial state and with the Bose factor 
in the final state, yielding: 



(27T) 

N 2 - 1 Tfila N 



■7- 9 2 T « 70, (9) 



AT g 2 h 2 4vr 

where we have made use of (|j). From this result, which has the same structure as the 
expression (|3|) for 7 , it is obvious that the plasmon mass m g as well as h cancel from the 
scattering rate. 

III. LYAPUNOV EXPONENTS 

The Lyapunov exponents measure the growth rate of infinitesimal perturbations around 
an exact solution of the classical lattice Yang-Mills equations. Since the maximal Lyapunov 
exponent Ao was shown to be independent of the lattice spacing, we assume that we can work 
in the continuum limit whenever adequate. If A^x, t) is an exact solution of the Yang-Mills 
equations, the linearized equation for a small perturbation a^(x,t) around is 

D 2 a^ - D^D u a u - 2i[F^ a v \ = 0. (10) 

Here D^(A) = — ilA^, } is the gauge covariant derivative where the bracket denotes the 
Lie algebra commutator, and is the field strength tensor associated with the background 
field Ap. 

The numerical approach to the determination of Ao proceeds by solving (|T0|) for an 
arbitrary initial condition a^x, 0) and measuring the growth rate of the norm of a^(x,t). 
To be precise, the maximal Lyapunov exponent was determined in [|T],0 from the logarithmic 
growth rate of the "distance" between neighboring field configurations, defined on the lattice 

as 

^M,U e ] = ^-J2\trU p -tvU± (11) 

where Ue are the group valued link variables, U p denotes the elementary plaquette operator, 
and N p is the total number of spatial plaquettes. In the continuum limit, the distance 



measure (11) takes the form 



V[A^ Afj\ oc J d 3 x |tr B'(x) 2 -trB 



(12) 



where B(B') are the magnetic fields associated with the gauge potential A^A'^). In going 
from (|TTD to (|12|) we have suppressed the constant factor (g 2 a/2N p ), since we are interested 
only in the growth rate of (lnT>). For an infinitesimal perturbation a M that is a solution of 
the linearized equation flTUD, we obtain: 



P[a M |A M ] = V[A^ + a^, A^} oc 



d x 



fdjti B 2 ) \ 1 ( d 2 (tr B 2 

tr l^AT a v + 2 tr IdAM; 



(13) 



The maximal Lyapunov exponent is then defined as 



hm Im.^o)^] 



(14) 



to— 1 P[a M (0)|A M ]' 
In practice, every randomly chosen initial configuration ^4^(0) with a fixed average energy 
density has been found to yield the same value for the maximal Lyapunov exponent Ao- The 
numerical calculations show that the maximal Lyapunov exponent depends only weakly 
on the lattice size and extrapolates smoothly to the limit of spatially homogeneous gauge 
potentials on a l 3 lattice. We take this as an indication that Ao is associated with long 
wavelength perturbations a M (x, t) in an appropriately chosen gauge. 



IV. ERGODIC LIMIT 

We now propose to make use of the fact, noted in the Introduction, that the background 
gauge field A^Xjt) rapidly approaches thermal configurations, by replacing the long-time 
average of the growth rate of (In T>) by the canonical average over background gauge fields 
A^, where the temperature T is chosen such that the thermal energy density equals the 
average energy density of the time-dependent background field A^x^t). The replacement 
of the temporal average by the canonical average relies on two conditions: The autocorre- 
lation function of the background field A^x, t) must decay on a time scale that is short 

7 



compared with the time to required for the calculation of the Lyapunov exponent, and the 
time evolution of the background field must be ergodic on the time scale t . 

The ergodicity of the background gauge field is assured by its dynamical chaoticity on 
time scales long compared to the inverse of the positive Lyapunov exponents, hence the 
second condition is fulfilled |16|| . On the other hand, if the first condition were violated, 
the Lyapunov exponent would depend on the starting configuration A^(x,t). In numerical 
studies [1-4] we have found that this is not the case. A direct study of the autocorrelation 
function performed by us has shown that the first condition is also satisfied. These conditions 
are in accordance with the g 2 T <C gT <C T hierarchy assumed in hot perturbative gauge 
theory. 

The maximal Lyapunov exponent is then obtained from the relation 



Ao~^ln(PMt)]) T , 

where the distance measure ( |i~3l) in a thermal background is 

'<9(tr B 2 )\ 



(15) 



(PM) T oc / d 6 x 



tr 



OA, 



(T) 



T 



1 

+ -tr 
2 



' <9 2 (tr B 2 ) \ 

, dApdA,, I rp 



(16) 



The first term in ( |T5| ) vanishes, because the thermal average of any quantity transforming 
under the adjoint representation is zero. In the second term, the thermal average projects 
on to the singlet part of <9 2 (tr B 2 )/dA fJi dA u , yielding 

'<9 2 (tr B 2 ^ 



«*H(S^) T tr( 



(17) 



Since the averaged value of T> is quadratic in the field fluctuations the Lyapunov 
exponent defined through the magnetic energy distance measure is twice as large as the one 
defined by the dominant exponential growth rate of the fluctuations of the elementary field 



A [^ = 2 t limilnJM^Jl 
to^°°t IK(o)|| 



(18) 
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V. CLASSICAL SPECTRAL FUNCTION 



Solving the classical equations of motion one deals with a problem essentially differ- 
ent from perturbative field theory: instead of investigating transition amplitudes between 
scattering states we follow the evolution of a given initial configuration from a time t — 
forwards. The appropriate method to analyze this evolution is not the Fourier transforma- 
tion as in quantum field theory, but the Laplace transformation. Its inverse transformation 
is then calculated along a path which has all poles of the spectral function on its same side; 
the path's position is shifted accordingly, compared with the Fourier transformation. 

The classical solution of the equations of motion for field perturbations therefore explores 
in forward time direction all poles of a free oscillator (or wave) equation. In case of chaotic 
Hamiltonian dynamics the solutions are both exponentially growing and damped giving rise 
to poles of the Laplace transform with positive as well as negative real parts. 

Making the formal connection between Laplace and Fourier transformation through a 
complex rotation of the frequency variable, s = icu, the inverse Laplace transformation path 
runs above all poles in the complex cu-plane. As a consequence in either case (oscillatory 
or chaotic) the integration path for the inverse Laplace transformation includes all poles for 
positive time and none for negative time while the Fourier transformation includes upper 
half plane poles for the advanced (negative time) and lower half plane poles for the retarded 
(positive time) propagator. 
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FIG. 1. The integration paths in the complex frequency plane for the inverse Laplace (L) and 
Fourier (F) transformations. 

The position of the poles obtained in a classical time — forward calculation may have 
in general both positive and negative imaginary parts. Therefore a better quantity for 
comparison is the spectral function which also considers poles in the whole complex em- 
plane. 

Summarizing this argument the position of all poles of a spectral function can be obtained 
from the linearized classical equations of motion for field perturbations (in the leading order 
of an h expansion), but the retarded and advanced propagators used to solve scattering 
problems in perturbative field theory discard the unsuitable ones due to their very definition. 
A positive Lyapunov exponent in Hamiltonian (energy conserving) dynamical systems, on 
the other hand, always occurs together with its negative counterpart — Liouville's theorem 
ensures it. Therefore studying positive exponential rates gives an information about the 
position of the poles of damped retarded and advanced propagators as well. 

The growth or damping rate, or the oscillation frequency of small amplitude fluctuations 
in a classical dynamical system is studied by linearizing the classical equations of motion. 
This procedure leads to a new differential operator whose spectrum gives the poles of the 
classical spectral function. Odd parity under time reflection, real valuedness and normaliza- 
tion conditions then determine the relative weights of the pole terms. 

The differential operator belonging to the linear perturbation propagation equation (10) 
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is identical with the second variation of the classical action, S"'[A], taken at the background 
field configuration A which is a solution of the classical equation of motion S'[A] = 0. 
Here the prime means variation with respect to A. Considering the generating functional 
of connected Green functions, the two point function is just the inverse of this differential 
operator, 

G[A,A'} = (AA') - (A)(A') = (S"[A})-\ (19) 

in the Gaussian approximation to the small amplitude fluctuations. So the linear perturba- 
tion propagation in classical equations of motion gives information about the saddle point 
approximated generating functional. 

Now aiming at the description of long wavelength plasmon damping we may neglect 
spatial derivatives and write the general form of the classical, linearized perturbation prop- 
agation equation (10) schematically as 

,; 2 



a(t) = 0. (20) 



The spectrum of this operator contains two poles on the real axis u = ±fi if Q 2 (t) is a positive 
constant. This case, familiar from zero-temperature perturbative field theory, describes small 
oscillations determining the real poles of the spectral function and the familiar retarded and 
advanced propagators. In classically chaotic, highly excited systems, however, it happens 
that tt 2 (t) is negative. This causes exponentially growing fluctuations — a typical source of 
chaotic behavior. 

In order to gain a qualitative understanding about the (classical) spectral function of 
chaotic systems we consider Q 2 (t) as a Gauss-distributed stochastic variable [|17|- It can 



have both negative and positive values, and its time variation is replaced by the ensemble 
variation due to the ergodic property of classically chaotic dynamical systems discussed in 
the previous section. In this limit the probability distribution of the frequency squares, 
P(Q 2 ), is determined by its two lowest moments, 
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(ft 4 ) - (ft 2 ) 2 = 4a 2 7 2 (21) 

parametrized by two real parameters a and 7. This parametrization reflects the fact that 
while (ft 2 ) can either be positive or negative, the width of its distribution is always positive. 

The stochastic average of the differential operator for the fluctuations has to be carried 
out on the quadratic level, because with the Gaussian distribution we assumed white noise 
property of the stochastic quantity. We get 

^(y 2 - ft 2 ) 2 ^ = u 4 - 2(ft V + (ft 4 ) = 
{uj — a — 27) (cu — a + i^){uj + a — 27) (cu + a + ij). (22) 

This result exhibits the symmetric four pole structure typical for a spectral function describ- 
ing classical plasma oscillations 

1/1 1 



Aid) 



Ainu \u — a — 27 u — a + 27 
1 1 
u + a + ry u + a — ry 



+ - -A-^) <*>) 



yielding the Lorentz shape 



A ( U ) = -7-5 o 27 ^9 ,,99 " ( 24 ) 

7T (u 2 - a 2 - 7 2 ) 2 + 47 2 u; 2 

The relative signs of the pole terms follow from the definition of the spectral function as 
the difference between the advanced and retarded propagators and from its odd time parity 

A(—u) = —A(u). The normalization factor 1 /2a; ensures that 

00 

J du uA(u) = 1, (25) 

—00 

so in each mode exactly one boson is counted by the spectral function A{uj) |L§] . 

This particular, four pole spectral function describes a general solution of the stochasti- 
cally averaged perturbation propagation equation which behaves like 

a(t) = Ae iat -^ + Be~ iat - Jt + Ce iat+Jt + De"^ 7 '. (26) 
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After some initial oscillations the exponential growth dominates the long time behavior of 
|a(t)|. It is exactly this, which has been seen in numerical calculations. The conclusion 
of this argument is that the Lyapunov exponent of elementary field fluctuations averaged 
ergodically is equal to the classical gluon damping rate as expressed by the imaginary part 
of the pole positions in the spectral function A(uS). 

We note that in a recent publication f|]I| a similar Gaussian model for the chaotic in- 
stability in general Hamiltonian flows has been investigated. Our result presented above 



recovers the more general one of ref |L9[ for vanishing expectation value of the noisy oscil- 
lator frequency square (a = 7) after substituting a characteristic timescale r = I/7 in the 



general formula (19) of ref. ||T9 



VI. CLASSICAL AND QUANTUM GLUON DAMPING 

Finally we argue again that the leading order gluon damping rate (3) obtained in hot 
perturbative QCD (pQCD) is classical, i.e. it retains its value in the classical limit h — > 0. 
This fact has been argued before in section II. Here we briefly reconstruct the argument 
and resolve some technical issues. This concludes our reasoning about the equality of the 
Lyapunov exponent of chaotic classical lattice gauge theory and the gluon damping rate at 
rest in a hot plasma. 

The gluon damping rate in hot pQCD is obtained from the definition (6) dividing the 
imaginary part of the self energy by the thermal gluon mass m g = gT/^/h. The general 
one-loop form of the self energy contains an integral over hard momenta, a factor of g 2 and 
the phase space distribution of thermal gluons 

ImllK, 0) = g 2 h J d A k ^n(k)f(k/m g ), (27) 

where the complicated algebraic expression f(k/m g ) depends only on scaled momentum 
variables. Using now the long wavelength approximation the phase space distribution of 
thermal gluons is replaced with its classical counterpart, n{k) rs T/hu, leading to 
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Imn(m S) 0) oc g 2 T J duf(u/m g ). (28) 

Scaling the integration variable with the Debye mass — which is of quantum origin contain- 
ing the Planck constant — we see that the imaginary part of the 1-loop gluon self energy in 
a hot plasma is proportional to m g . It follows that the gluon damping rate obtained using 
"classical" thermal gluons does not depend on the Debye mass and Planck's constant, 

j = g 2 Tj dxf(x), (29) 

showing that the result (3) is essentially classical. 

Finally, it is still to show whether non-pole contributions to the self energy in the field 
theoretical calculation do not interfere with the above arguments. The one-loop spectral 
function used there as an input contains a pole term picking up the zeroes u(k) of the 
inverse propagator corresponding to collective plasma modes to the lowest order and a cut 
term describing the effect of scattering on thermally excited spacelike modes: 

A(k, u) = Z{k)5{uj 2 - u{k) 2 ) + p(k, u)0{k 2 - u 2 ). (30) 

The cut coefficient j3(k,u) is related to the real and imaginary parts of the self energy 
n(Jfe,w): 

±imn , 

0(k,u) = E 5 ~. 31 

V ; (P-cu 2 + Ren) 2 + (ImIT) 2 V ; 

The respective self energies for the transverse and longitudinal excitations to leading order 



in hot perturbative QCD are [20 



U t (k, u) = m 2 x 2 ( 1 + ^^!hi^^ + -m 2 x(l - x 2 ) (32) 

and 

U e {k,u) = k 2 + m 2 (2 -x ln^^ ) - mm 2 x (33) 

with x = io/k and m 2 = 3m 2 /2. Using these forms one obtains the following cut parts of 
the retarded Fourier transform of the spectral function, A(t, k), for small k/m 
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A cut , t (M) -> -26(t)-^ 

rrr 

1 , 4x(l — x 2 ) sin ktx , . 

cix i '— (34) 

4x 4 (l + i^lniif) + x 2 (l-x 2 ) 2 



A cuM 0, t) -> -26(t)— -J /' rfx xsmktx _ (35) 



and 



2-xln^ +x 

Since the integrand is bounded, the cut contribution cannot grow exponentially with time 
and hence does not contribute to the maximal Lyapunov exponent (|i~5|). In fact, the cut 
contribution vanishes in the long wavelength limit k — > 0. 

This leaves us with the pole part, which remains finite in this limit. 



VII. SUMMARY 

This concludes our argument establishing a connection between the classical Lyapunov 
exponent and the gluon damping rate in hot perturbative QCD. We note that some elements 
of the argument are heuristic, in particular, the replacement of the long-time average of the 
growth rate of fluctuations around a specific field configuration by the thermal average. This 
reasoning assumes that the growth rate, or equivalently the plasmon damping rate, depends 
only on coarse-grained properties of the gauge field. We believe that this is so, because the 
one-loop calculation of the damping rate 70 only involves soft loop momenta || and hence 
does not depend on details of the short- distance fluctuations of the gauge field. 

Because of the general nature of our argument, we conjecture that the complete spec- 
trum of Lyapunov exponents obtained in || reflects the spectrum of damping rates j(k) 
of excitations in a thermal bath. If this were true, it would confirm our assumption that 
j(k) < 7o- Since, at present, it is not known whether "f(k) is a quantity with a classical 
limit for k 7^ 0, the identification with the Lyapunov spectrum remains a conjecture. We 
finally note that if the correspondence between ergodic and canonical averages holds up for 
other physical quantities, transport coefficients of nonabelian gauge fields at the classical 

15 



scale {g 2 T), such as magnetic screening 21] or color diffusion [22], could possibly also be 
calculated by real-time evolution of classical gauge fields on a lattice. 



Acknowledgements: We thank U. Heinz, S. G. Matinyan, H. B. Nielsen, G. K. Savvidy, 
and M. Thoma for illuminating discussions and the referee for pointing out a mistake in 
our original manuscript. This work was supported in part by the U.S. Department of En- 
ergy (grant de-fg05-90er40592) and in part by the Collaboration Agreement between the 
Norwegian Research Council (NFR) and the Hungarian Academy of Science (MTA) (grant 
422.92/001). One of us (B.M.) thanks the Physics Department of Tokyo Metropolitan Univer- 
sity, especially H. Minakata, for their hospitality and support during his visit there. T.S.B. 
acknowledges the hospitality and the support of Physics Department of Bergen University 
and of L. P. Csernai during his visit there. 



16 



REFERENCES 

[1] B. Miiller and A. Trayanov, Phys. Rev. Lett. 68, 3387 (1992). 
[2] C. Gong, Phys. Lett. B298, 257 (1993). 
[3] C. Gong, Phys. Rev. D49, 2642 (1994). 

[4] T. S. Biro, C. Gong, B. Miiller, and A. Trayanov, Int. J. Mod. Phys. C5, 113 (1994). 
[5] E. Braaten and R. D. Pisarski, Phys. Rev. D42, 2156 (1990). 

[6] V. V. Klimov, Zh. Eksp. Teor. Fiz. 82, 336 (1982) [Sov. Phys. JETP 55, 199 (1982)]. 
[7] H. A. Weldon, Phys. Rev. D26, 1394 (1982). 

[8] E. Braaten and R. D. Pisarski, Nucl. Phys. B337, 569 (1990); Phys. Rev. Lett. 64, 
1338 (1990). 

[9] J. Frenkel and J. C. Taylor, Nucl. Phys. B334, 199 (1990). 

[10] J. P. Blaizot and E. Iancu, Nucl. Phys. B390, 589 (1993); Phys. Rev. Lett. 70, 3376 
(1993); Nucl. Phys. B417, 608 (1994). 

[11] U. Heinz, Ann. Phys. (NY) 161, 48 (1985); ibid 168, 148 (1986). 

[12] J. C. Taylor and S. Wong, Nucl. Phys. B346, 115 (1990). 

[13] R. Efraty and V. P. Nair, Phys. Rev. Lett. 68, 2891 (1992); Phys. Rev. D47, 5601 
(1993). 

[14] The sign conventions in ref. 5 are such that a positive imaginary part of *n implies a 
negative imaginary part of the pole energy. In writing (6) we have already continued 
*n t to Minkowski space. 

[15] P. F. Kelly, O. Liu, C. Lucchesi, and C. Manuel, Phys. Rev. Lett. 72, 3461 (1994). 

[16] Note that the infrared limit of the Yang-Mills equation has been shown not to be 

17 



completely ergodic, with some very small regions of stable trajectories remaining in 
phase space. See: P. Dahlquist and G. Russberg, Phys. Rev. Lett. 65, 2837 (1990); 
R. Marcinek, E. Pollak, and J. Zakrzewski, Phys. Lett. B327, 67 (1994). It is unknown, 
if such a phenomenon persists in the case of spatially dependent Yang-Mills fields. 

[17] see e.g. N. G. van Kampen, Phys. Rep. 24, 171 (1976). 

[18] G. Baym and N. D. Mermin, J. Math. Phys. 2, 232 (1961). 

[19] L. Casetti, R. Livi and M. Pettini, Phys. Rev. Lett. 74, 375 (1995). 

[20] R. D. Pisarski, Physica A158, 146 (1989). 

[21] A. D. Linde, Phys. Lett. B96, 389 (1980); T. S. Biro and B. Muller, Nucl. Phys. A561, 
477 (1993). 

[22] A. V. Selikhov and M. Gyulassy, Phys. Lett. B316, 373 (1993); H. Heiselberg, Phys. 
Rev. Lett. 72, 3013 (1994). 



18 



